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Abstract 

We compute NNLO {0{al)) corrections to the hard-scattering kernels entering the QCD 
factorization formula for B V'y decays, where V is a light vector meson. We give 
complete NNLO results for the dipole operators and Qg, and partial results for Qi valid 
in the large-/?o limit and neglecting the NNLO correction from hard spectator scattering. 
Large perturbative logarithms in the hard-scattering kernels are identified and resummed 
using soft-collinear effective theory. We use our results to estimate the branching fractions 
for B —y K*^ and Bg —>■ 07 decays at NNLO and compare them with the current 
experimental data. 
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1 Introduction 



Radiative B V'~) decays, where is a light vector meson, are processes of particular 
interest in flavor physics which are already accessible at the 5-meson factories at SLAG 
and KEK; current measurements [1-7] yield the branching fractions presented in Table [H 
These decays provide independent constraints on the shape of the unitarity triangle, 
determining the side Rt of this triangle through the ratio of branching fractions for B 
{p,uj)^ and B — »• K*^ decays. This information is complementary to the constraints on 
the ratio of CKM matrix elements |V^d/V^s| obtained from the recent CDF measurement 
of the mass difference AM,, in the Bg — Bg system [8] and the already precise knowledge 
of the Bd — Bd mass difference AMd [7]. Moreover, measurements of the CP-asymmetries 
in B ^ (p, uj)'~f decays and the isospin-violating ratio of the charged and neutral B p'y 
modes would determine the inner angle a of the unitarity triangle. 

The calculation of the branching fractions for B V'y decays requires the evaluation 
of the hadronic matrix elements of the operators in the effective weak Hamiltonian. For 
B decays the weak Hamiltonian is [9, 10]: 



p=u,c 



CiQi + C2Q2 + CiQi 



1=3 



where A^,''^ = V^^Vpb (unitarity of the CKM matrix implies that aJ^^ = -(a1^^ + Ai''^) and so 
contributions from diagrams with top-quark loops are included implicitly). The relevant 
four-quark operators Qi and Q2 are 

Qi = {qp)v^a ipb)v-A, Ql = {qiPj)v^A iPjbi)v-A, (2) 

and the electromagnetic and chromomagnetic penguin operators and Qs are 

Q7=-^^gcr-ni + 75]&F,., Qs = -^-^qa^'[l + l,]T^bG;^. (3) 

Here q = d or s, and the convention for the sign of the couplings corresponds to the 
covariant derivative iD^^ = idfj, + eQjA^ + gT^-A^j^, with and representing the 
photon and gluon fields respectively, and Qe = —1 etc. The factor fnb{fi) is the MS mass 
of the b quark. The Wilson coefficients Ci have been known within the next-to-leading 
logarithmic approximation (NLL) for over a decade (for a review, see [11]), and have 
been recently calculated at next-to-next-to-leading logarithmic order (NNLL) in a series 
of papers [12-16]. In the present work we focus on the most phenomenologically relevant 
operators, which are Qi, Qi, and Qs- The matrix elements of the QCD-penguin operators 
Qs, . . .Qq first contribute at 0{as) and are multiplied by small Wilson coefficients in the 
weak Hamiltonian Tieff (III)- The contribution from Q2 starts at (9(q;^). 

It has been shown that in the heavy-quark limit a factorization framework (called QCD 
factorization) can be applied to B V'-f decays [17-23] (see [24,25] for phenomenological 
updates to NLO, and [26,27] for the alternative "perturbative QCD" approach). In 
particular, the matrix element of a given operator in the effective weak Hamiltonian can 
be written in the form 

{V^ \Qi\ B) = F^-^^^ Tl + Jdoo du <Pl{u) T^{uj, u) . (4) 
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Table 1: Status of the fi-meson radiative branching fractions (in units of 10 ^) from the 
BABAR, BELLE and CLEO collaborations and their averages by HFAG [7]. The entry 



for 



is from the recent BELLE measurement [6]. 



Mode 



BABAR 



BELLE 



CLEO 



37.6+1^ ± 2. 



HFAG 



38.7 ±2.8 ±2.6 
39.2 ±2.0 ±2.4 



42.5 ±3.1 ±2.4 
40.1 ±2.1 ±1.7 



45.5 



+7.2 



±3.4 



40.3 ±2.6 
40.1 ±2.0 



:oil ± 0.06 
o.4o;g:^g ± 0.05 



n i;c;+U.42+U.09 
'-'■"^^-0.36-0.08 
1 or +0.37+0. 07 
-■-■^^-0.33-0.06 



0.93l°i^ 



B+ 
B°- 

50- 



> 



0.79^ 



0.56 



+0.34+0.05 
-0.27-0.10 



< 13 

< 17 

< 9.2 



B 
B 



40.4 ±2.5 
1.25 ±0.25 ±0.09 



42.8 ±2.4 
1.32 



+0.34+0.10 



-0.31 
TTST 



-0.09 
T2 



43.3 ±6.2 
< 14 



41.8 ± 1.7 
1 90+0.31 

-'-•^°-0.29 



5., 



07 



57 



-15-17 



The non-perturbative effects are contained in F^^^-^ , the B ^ V transition form factor at 
= 0, and in (j)^ and 0^, the leading-twist light-cone distribution amplitudes (LCD As) of 
the B- and V^-mesons. The hard-scattering kernels T/ and T/^ include only short-distance 
effects and are calculable in perturbation theory. Contributions to the kernel are closely 
related to the virtual corrections to the inclusive decay rate, and are referred to as vertex 
corrections. Those to the kernel are related to parton exchange with the light quark 
in the S-meson, a mechanism commonly referred to as hard spectator scattering. It is 
expected that the factorization formula is valid up to corrections of (9(Aqcd/''^6)- 

The derivation of the factorization formula from a two-step matching procedure in soft- 
coUinear effective theory (SCET) [28-31] has provided additional insight into its structure. 
An advantage of the effective field-theory approach is that it allows for an unambiguous 
separation of scales and an operator definition of each object in the factorization formula. 
The technical details for B V'~f have been provided in [32-34]. In the SCET approach 
the factorization formula is written as 

{Vj IQilB)^ AiC^Cv^ + ^I^eIIx^ J du du 4>l{u) 4>l{u) tY{u, u) , (5) 

where F and fv^ are meson decay constants. The SCET form factor Cv^ is related to 
the QCD form factor through pcrturbativc and power corrections [35,37-42]. In SCET 
the perturbative hard-scattering kernels are the matching coefficients AjC"^ and . They 
are known completely to next-to-leading order (NLO) {0{as)) in renormalization-group 
(RG) improved perturbation theory [34]. In this paper we make steps towards a complete 
analysis at next-to-next-to-leading order (NNLO) by obtaining full results for the hard- 
scattering kernels for the dipole operators Q-j and (^8, and partial results for Qi, valid in 
the large-/3o limit and neglecting NNLO corrections from spectator scattering. 

The hard-scattering kernels are found by matching certain partonic matrix elements 
in QCD with those in the effective theory. For the vertex corrections the relevant matrix 
elements are {s^\Qi\h). The loop corrections in the effective theory can be made to 
vanish by matching on-shell, so the main obstacle is the evaluation of the QCD matrix 
elements. However, these matrix elements are just the virtual corrections to the inclusive 
B — >■ Xg^ decay rate. Exact results to 0{a1) were obtained for Q-j in [43,44] and for 
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in [45]. For Qi the virtual corrections at 0{as) were calculated in [46-48], but those at 
0{al) are known only in the large-/3o limit [49,500. A calculation that goes beyond this 
approximation by employing an interpolation in the charm quark mass uic was reported in 
[52] , and has been used in estimating the NNLO branching fraction for the inclusive decay 
B — »• [53] . However, as the calculation was not split into virtual and bremsstrahlung 
contributions, those results cannot be used in the SCET matching calculation. Therefore, 
while we obtain exact NNLO results Q-j and Qs, for Qi we are restricted to the large- /5o 
limit. Our results provide an explicit check on factorization at NNLO. 

Corrections from spectator scattering are included in the hard-scattering kernel t^^ and 
first contribute to the branching fraction at NLO {0{as))- A complication of spectator 
scattering is the presence of two widely separated perturbative scales ml S> mftAgcD- The 
SCET approach provides a systematic framework for separating contributions from these 
two scales. In SCET the hard-scattering kernel for a given operator is sub-factorized 
into the convolution of a hard-coefficient function with a universal jet function, in the 
form ^ 

tl\u,u;)= [ drA,C^i(r)j^(r,n,^) = A,C^i^j^. (6) 
Jo 

The hard coefficients AiC^^ contain physics at the hard scale rrih, while the jet function 
j± contains physics at the hard-collinear scale ^/mbA. The hard coefficient is identified 
in a first step of matching QCD — > SCETi, and the jet function in a second step of 
matching SCETj —>■ SCETn. Details have been worked out for B — »• V'j in [32,34], for 
heavy-to-light form factors in [36-42], and for B — > PP in [54,55]. 

The effective field-theory techniques are crucial for providing a field-theoretical defi- 
nition of the objects in i^, and for resumming large perturbative logarithms of the ratio 
mb/AqcD in the t]^. In the effective-theory approach resummation is carried out by solv- 
ing the renormalization-group equations for the matching coefficients AjC^^. Since these 
coefficients enter the factorization formula in a convolution with the jet function j±, their 
anomalous dimension is a distribution in the variables r and u. The evolution equations 
must be solved before performing the convolution with j^. Therefore, resummation is not 
possible in the original QCD factorization formula (jlj), where the hard-scattering kernels 
T/^ are obtained only after this convolution has been carried out. 

While the SCET formalism is indispensable for resummation, in the actual matching 
calculations one can also use the diagrammatic method of expanding by regions [56] in 
order to separate hard from hard-collinear effects as in ([H]). This method was used to 
analyze loop corrections to spectator scattering for the case of the B ^ it form factor 
in [57], and for B vrvr in [58]. In both cases the results were shown to be equivalent to 
those obtained directly in SCET. We use similar techniques here to compute the NNLO 
correction from the hard-scattering kernel t^g. Our result for the one- loop correction 
at the hard-collinear scale agrees with explicitly confirming the universality of the 
jet function predicted by SCET. Since the NNLO corrections from are known from 
the form factor analysis [40,42], the main obstacle to a complete treatment of spectator 
scattering is the NNLO matching calculation for Qi. 

The paper is organized as follows. In Section [2] we explain the SCET factorization 
framework and define the hard-scattering kernels. The SCET matching calculations are 

"^these results are obtained by calculating the 0{alnf) terms and then replacing Uf — 3/3o/2, 
according to the hypothesis of "naive non-abelianization" [51]. 
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carried out for the vertex corrections in Section [3] and for the hard-spectator corrections 
in Section HI In Section [5] we describe the numerical analysis and estimate the branching 
fractions for B —>■ K*'~f and Bs <f)'~f decays at NNLO, comparing our results with the 
current data and identifying the theoretical uncertainties. We conclude in Section [HI 
Results for the partonic matrix elements taken from calculations for inclusive B —>■ X^'-f 
decays are relegated to the Appendix, along with some of the SCET matching functions 
obtained in previous work and details of the renormalization-group analysis. 



2 Factorization and the hard-scattering kernels 

In this section we explain the closely related issues of factorization and extraction of the 
hard-scattering kernels. The objects of interest are the hadronic matrix elements 

An analysis in [34] used a two-step matching procedure in SCET to show that these 
matrix elements can be written in the form ([5]) to all orders in perturbation theory and to 
leading order in 1/m;,. In this paper we work out a large set of effective-theory matching 
coefficients at NNLO in perturbation theory. These are obtained by replacing the hadronic 
states by partonic ones and calculating the matrix elements in perturbative QCD. Showing 
that the partonic rate can be brought into the form ([5]) demonstrates factorization and 
provides expressions for the hard-scattering kernels. 

To calculate the partonic matrix elements requires the evaluation of multi-scale Feyn- 
man integrals. It is advantageous to perform these integrals using the method of re- 
gions [56]. This not only provides a simple way to obtain results at leading order in l/rrih, 
but also a factorization of momentum scales at the level of Feynman diagrams. In this 
method the loop integrations are split into a sum of different regions, in which the loop 
momenta satisfy a fixed scaling. This allows for a Taylor expansion under the integral in 
each region, which is subsequently integrated over all space. The integrals are performed 
in dimensional regularization, where scaleless integrals are set to zero. The sum of the 
results for all the regions recovers the full integral, expanded in l/rub. 

A number of different momentum regions appear in the analysis, both perturbative and 
non-perturbative. To identify these we first introduce two light-like vectors n± satisfying 
n^n_ = 2. We choose the outgoing vector meson to travel along the n_ direction, and 
define such that the velocity of the b quark is given by 

"2+2 ^ ' 

This definition implies v± = 0, and we shall always work in the reference frame where 
n_v = n^v = 1. To perform the expansion in l/rrif,, we define the parameter A^ = 
{Pb — rribvy and the dimensionless parameter A = A/rrih <C 1. The regions are classified 
according to the scaling of their light-cone components with the expansion parameter 
A. Denoting the light-cone components of a generic four-vector p by {n+p,p±,n-p), the 
relevant momentum regions are [34]: 
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Perturbative 



hard Tnb{l, 1, 1) 

hard-collinear /«f,(l, VX, A) 

Non-pert urbative 



soft ^b(A, A, A) 

collinear mfe(l,A,A^) 
soft-coUinear mf,(A, A^/^, A^) 

The connection between the SCET analysis and perturbative QCD is provided by the 
method of regions. In the effective theory, contributions from the perturbative regions 
are encoded in Wilson coefficients of operators built from fields representing the regions 
of lower virtuality. It is convenient to factorize the two perturbative scales ml and mfoA 
using a two-step matching procedure QCD SCETj SCETn- 

In the first matching step the hard scale ml is integrated out by matching the operators 
Qi onto a set of operators in SCETj. The effective theory SCETj involves fields for the 
hard-collinear and non-perturbative modes, multiplied by Wilson coefficients related to 
the hard region. For the case of B ^ V'j, the matching takes the form [34] 

Q, ^ A,C^J^ + A,C^i ^ + A,C^2 ^ J^l (8) 

The -k denotes a convolution over momentum fractions, as in ([6]). The momentum-space 
Wilson coefficients depend only on quantities at the hard scale m^. The exact form of the 
operators J*-*^ along with the relevant SCET conventions can be found in [34]: 

= {^W,,)^^{l-j,)K, (9) 

J""' = (eW^hc)^±AeJl+75)/i., (10) 

J""^ = {^Whc)AcJ±ii + i5)K- (11) 

Here e± is the polarization vector of the on-shell photon. The operators contain a hard- 
collinear quark field ^, a composite object Ahc, which in light-cone gauge is the hard- 
collinear gluon field, and Whc, a Wilson line. In SCET the 6-quark field is treated as in 
HQET. We have suppressed the arguments of the fields above, but must keep in mind that 
due to the non- locality of SCET the objects {^Whc) and ^hcj_ are evaluated at different 
points along the n+ light-cone, whereas is multipole expanded and evaluated at a point 
on the n_ light-cone (see, e.g., [31]). The -B-type operators are actually power suppressed 
in SCETi, but contribute at the same order as the A-type operator upon the transition 
to SCETn [36-38]. 

The matrix element of the operator is proportional to the SCET form factor (v±- 
The Wilson coefficients AiC^ multiplying this matrix element can be extracted from 
calculations in the inclusive B — > Xgj decay. Details are given in Section [3l In contrast 
to the QCD form factor, the SCET form factor contains no piece which can be written 
in the form of a (convergent) convolution of a hard-scattering kernel with the meson 
LCDAs [37, 380. The relation between the QCD form factor and the SCET form factor 

^although see [59] for a renewed discussion of this point. 
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is determined by the factorization formula [35,37,38] 



pB^v, = C^Jy^ + ^^f^- I rfu;rf«0»05:(«)t^i(c^,«). (12) 

Since the matrix element of is proportional to the form factor, the coefficient functions 
Cy^ and ty^ at NNLO can be determined from the results for Qj. The exact relation is 
given in (IHOl) below. 

The operators J^^^^ can be further matched onto four-quark operators in SCETn- For 
B V'j decays, only the operator J^^ is relevant. The matrix element of the four-quark 
operator onto which it matches factorizes into a product of LCDAs for the B and V 
mesons. The operator J^^, on the other hand, matches onto a four-quark operator whose 
renormalized matrix element has no projection on the pseudoscalar i?-meson LCDA. In 
matching the operator J^^ onto SCETn the hard-coUinear scale m^A is integrated out, 
and the associated Wilson coefficient is the jet function The final low-energy theory 
SCETn contains only soft, collinear, and soft-collinear fields. Factorization means that 
soft fields are restricted to the 5-meson LCDA, and collinear ones to the V^-meson LCDA. 
Since these two pieces communicate only through soft-collinear interactions, factorization 
amounts to showing that such contributions decouple from the hadronic matrix element of 
the SCETii operator. This was done in [34] . Thus the matrix element of the operator onto 
which J^^ matches is exactly of the form of the second piece of ([5]), with tf = AiC^^ -kj^. 
This same jet function appears in the factorization formula ( IT2l) for the form factor, where 
ty^ = Cy^ * J±- We can summarize this discussion by the following factorization formula 



This formula relates the hard-scattering kernels AiC^ and tf in to the Wilson coef- 
ficients from the two-step matching procedure in SCET, and provides a connection with 
the original formulation (jlj). For instance, using that A7C* ~ Cy , one can verify that 
Qj contributes to both terms in the SCET formulation, but only to the vertex term in 
the original formulation. 

A main result of our paper is an expression for the 0{af) correction to the hard 
coefficient AgC^^. We obtain it with a straightforward diagrammatic analysis using the 
method of regions, without the explicit formulation of SCET or the use of its Feynman 
rules. Since contributions from J^^ can be uniquely identified by the Dirac structure of 
the four-quark operator onto which it matches, the sub-factorization of the hard-scattering 
kernel into a convolution of a jet and hard function can be performed by separating out the 
contributions of the hard and hard-collinear regions multiplying this structure. Details 
are given in Section HI 

3 Vertex corrections 

We begin with the vertex corrections, extracting the contributions of the operators Qi, Qj, 
and Qs to the SCET Wilson coefficient at NNLO To do so we calculate the 
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partonic matrix elements 

{Q,) ^ {q{p)l{q)\Qi\h{p,)) 

to this same order in both SCET and QCD. This matrix element is chosen because it 
contains no external gluons and so matches directly onto the operator in ([H]). The cal- 
culation is performed with on-shell external quark states and both UV and IR divergences 
are regularized dimensionally. In that case the matching calculation is simple, because the 
loop corrections in SCET are scaleless and vanish. The matrix element of J"^ is just the 
tree expression plus counterterms from wave-function and current renormalization. The 
QCD matrix elements can be read off from the virtual corrections to the inclusive decay 
B Xg'^. Using that the on-shell wave-function renormalization factors in the effective 
theory are unity, and replacing the bare SCET current by its renormalized one, we have 

{Q^) = A(Q7,trcc) = A.C^Zj ( J.^J . (14) 

Here the Di are the scalar amplitudes in QCD, the AjC^ are the contributions of a given 
operator to the SCET matching coefficient, and Zj is the renormalization factor of the 
SCET current operator J^. Each of these quantities is determined as a series in a^- For 
the operators Qy^g we can obtain complete results at NNLO, while for Qi we can only 
provide an estimate using the large-/3o limit. 

We first consider tree level, where only Q-j contributes. For on-shell matching the 
spinors in QCD and SCET are equal to one another and we find 

471-2 ' ^ ' 

where the photon energy is 2E^ = 771^(1 — my/m%) ~ rub in the heavy-quark limit. At 
higher orders the matching coefficients can be read off from the functions according 
to the relation 

A,C^(mb, /i) = A7C^(°) lim Zj\e, m,, /i) A(e, m^, /x) . (16) 

The SCET current renormalization factor Zj is determined by requiring that the Wilson 
coefficient be free of IR poles. 

Before giving results for the higher-order corrections, we pause to explain a subtlety in 
the matching which first appears at two loops. The on-shell matrix elements of the QCD 
operators Qi are calculated in MS renormalization in the five- flavor theory, nf = ni + Uh 
with Uh = 1 for the b quark. However, in SCET 6-quark loops are absent and the matrix 
elements are calculated as an expansion in the four-flavor theory. In order to perform 
a correct matching, it is necessary to express the UV renormalized results in the five- 
flavor theory in terms of the four- flavor parameters of SCET. A similar problem arises 
when integrating out the top quark to match the Standard Model onto the effective weak 
Hamiltonian. The solution is to renormalize the coupling constant in the Uf = rih + ni 
flavor theory according to a^^^^ = Z'^''~^"''as, with (see e.g. [13,60]) 



112 2 
^CA--nf + -nh{l-N, 



(17) 



The function A^^ is fixed such that as is the MS-renormalized coupling in the four flavor 
theory. Its value is 

^(e)=e^^(^)V(l + e). (18) 
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Results for the scalar amplitudes Di in this renormalization scheme can be obtained from 
the MS results given in the Appendix by making the replacement 



as 4 
' 47r 3 



TT 



L + e L' + — 



12 6 



+ . . . 



where L = In fi/rrih. Note that this is just the standard decoupling relation when evaluated 
in four dimensions. 

We now give results for the Wilson coefficients, which we write in the form 



Air 



47r 



(20) 



We begin with Qj. Results can be given analytically, but since those for are only 
known numerically we treat Q-j the same. Using the scalar functions given in the 
Appendix we find 

AyC^^i) = Cf [-2L^ - 5L - 2Lqcd - 6.8225] , 

A7C^(2) ^ C2 (2L^ + 14L3 + 38.1449^2 + 56. 14711L + 7.8159) 

+CfCa (-4.8889L3 - 33.9758^^ - 92.3415L - 83.8866) 

+CFni (0.8889L^ + 6.8889^^ + 19.9050L + 23.8254) 

+CFnh (-1.3333L2 + 2.8889L - 0.810288) , (21) 

where one is to use = 4 and n/j = 1 in the above equation. In the one-loop result we have 
distinguished the logarithms Lqcd = In /iqco/'^fe and L = In fi/mf,. The /iqcd dependence 
cancels against the scale dependence in the effective weak Hamiltonian, whereas the n 
dependence cancels against the scale dependence of the SCET soft function (vi_ and the 
running coupling constant. At one loop it is straightforward to separate the logarithms 
by identifying the UV and IR poles in the individual Feynman diagrams. At two loops 
the distinction can be made by using the renormalization-group equation fl^^ below. We 
give explicit results for the case where L is distinguished from Lqcd in the Appendix, but 
in this section we quote the NNLO results only for L = Lqcb- 

We can use our results to determine the anomalous dimension of the operator up 
to two loops. The anomalous dimension is obtained from the coefficient zf^ of the 1/e 
pole term in the current renormalization factor and has the form 
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da.. 



Zf\mb,fi) = -rcusp(as) ln— + Y{as) 

rub 



(22) 



where Fcusp is the cusp anomalous dimension appearing in the renormalization-group 
theory of Wilson lines [61] (it has recently been calculated to three loops [62]; the result 
is listed in the appendix). The result for the renormalization factor to two loops is 

2L' 



+Cf - 



An, 



CfOs 


1 


5 


An 




2e 






0.5Cir 




r 









- - 2.5Cf + 2.75Ca - 0.5ni - 2CfL 



- 3.125Cf + 3.5447C^ - 0.5556/1; - 2CfL^ + {-5Cf + 3.6667C^ - 0.6667ni)L 



+- - 2.6525Cf - 3.4386Ca + 1.9799/1; + 



-4.1546C^ + l.lllln;)Lj 



(23) 
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from which we find 

7^ = ^(-4^-5) (24) 

+Cf 7^ ((-16.6183C^ + AAAAni)L - 10.6102Cf - 13.7545C^ + 7.9195^;) . 
\47r/ 

This is consistent with ([22]) and defines 7*^. The one- loop result was first obtained in [29]. 
We note that in this case the dependence in the renormalization factor Zj drops out 
after using ( |T9i) . This must be the case, since in the effective-theory current the h quark 
is integrated out and so its anomalous dimension cannot depend on Uh- Our result for 
the anomalous dimension, along with the relation 

/i-^A,C^ = 7^A,C^, (25) 

allows us to perform the separation of UV and SCET logs in the Wilson coefficients given 
in the Appendix. 

The same SCET current also appears in the study of the inclusive B —>■ Xg'j decay 
spectrum with a cut on the photon energy [63]. A result equivalent to our two- loop 
matching coefficient AyC"^^^-* with = was recently obtained in [64]. Translating our 
expression into the two-loop result for h{l) given in [64], we find numerical agreement. 
The dependence on n/j not taken into account in that work is negligible numerically. We 
can also check the two-loop anomalous dimension by using RG-invariance of the inclusive 
decay rate along with the anomalous dimensions of the jet and soft functions calculated 
in [65,66]. Here again the results agree. 

We repeat the calculation for Qg. In this case the one- loop result is IR finite. The 
two-loop matching equation also becomes IR finite after the results are expressed in terms 
of the renormalized current calculated above. This is a check on the effective-theory 
construction, according to which the IR poles in the QCD amplitudes for each operator 
in the weak Hamiltonian are absorbed by the same SCET current. For the coefficient 
functions we find 

AgC^W = [2.6667Lqcd + 1-4734 + 2.0944i] , 
AgC^^^^ = [5.3333L3 + 32.2802^2 + 50.9612L+ 1.8875 

+i(4.1888L2 + 31.4159L + 29.8299)] 

+CfCa [15.1111L=^ + 31.6617L + 2.3833 + i(23.7365L + 28.0745)] 

-Cpni [1.77781.2 + 4.0386L + 1.7170 + i(2.7925L + 4.4215)] 

+CFnh [1.7778L2 - 2.0741L + 0.8829] . (26) 

Finally, we consider the four-quark operators Qi and Q2- At NLO the contribution 
from Qi can be obtained as an expansion in ml/m^, whereas that from Q2 vanishes. 
To extract the NNLO results for these operators would require the QCD amplitudes Di 
and D2 to this same order, which involves the calculation of a large set of three-loop 
graphs. These corrections are known exactly only in the large-/5o limit, in an expansion 
in z = m'^/ml [49]. Within this approximation the result for Q2 vanishes, and that for 
Qi can be written as 

A.C^W = [-3.8519LQCD + r^'\z)] , (27) 

AiC^(2) = -^^Cf \2AQ91L' + l^'\z)L + r^'\z)] , 
2 nif, 
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where we have replaced nf —>■ — 3/5o/2 as appropriate in the large-/5o hmit. Within this 
hmit it is also consistent to set the ratio rrib/rnb to unity, as we shall do in the numerical 
analysis of Section [5l Since in the large-/3o limit the amplitude is IR finite, we can read 
off the functions r*^*^ and Z*^^^ directly from the results for inclusive B — > X^'j decay. 
Converting to our notation we have 

r^'' = ^, r(^) = ^ I''' = -jr. (28) 

('2'! ('2) 

where r2 is defined in eq. (2.35) of [46], and r2 , ^2 sq. (22) of [49]. As an example, 
for mcjvrib = 1.2/4.8 we have 

^^^A{i) ^ ^(j^ [-3.8519LQCD - 3.4529 - 0.5138i] , (29) 
niij 

AiC^(2) ^ [2.4691^2 + 4.9083L + 5.1203 + z(0.9953L + 1.6014)1 . 

2 nib 

There are two major uncertainties associated the large-/3o limit. The first is that there 
is no way to quantify the size of the terms in AiC"^'-^^ not captured within this limit. The 
second is that the higher-order calculation does not resolve the perturbative ambiguities 
in the ratios of quark masses mb/frib and nic/rnb in the lower-order coefficient AiC^^-^^: 
the difference between mass renormalization schemes in these ratios is a correction pro- 
portional to Cpcts and set to zero in the large-/3o limit. We discuss these uncertainties in 
more detail in the numerical analysis of Section O 

In Section [3] we will be interested in the dependence of the branching fractions on 
the choice of renormalization scales. Both the matching coefficients AiC^ and the SCET 
soft function ^vi depend on the SCET factorization scale fi. It is convenient to use 
the renormalization group to determine the coefficients AiC^ at an arbitrary scale fi, 
given their value at a matching scale fih ~ /^qcd ~ ^b- This allows us to fix /i = nib 
and determine the soft function (vj_ only at this single scale. We can then study the 
dependence of the branching fractions under variations in fih and /iqcD, under which it 
is formally invariant. The relevant RG formalism was worked out in [41]. The expression 
we need is 

AiC^{mb, /ife, /i) = — exp[S{nh, /i) + ajinh, n)]AiC^{mb, Hqcd = l^'h, l^'h) ■ 

(30) 

In the above equation we have correlated the scales /xqcd = /^h for simplicity, although 
we can keep them separate using the results in the Appendix. With this choice the 
dependence on fih = A'-qcd on the left-hand side cancels against the dependence in the 
effective weak Hamiltonian, so that the branching fractions are invariant under variations 
of the matching scale fit- The RG exponents S and a, and aj are given by 

a(/ii,/i2) = / -r^rcusp(a), (32) 
aj(/ii,/^2) = / -T-^7j(a)- (33) 
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These exact solutions are evaluated by expanding the anomalous dimensions and the QCD 
/9-function as perturbative series in the strong coupling. We can do this to two-loop order 
for a and aj, and to three-loop order for 5*. The expansions to this order are listed in the 
Appendix. 



4 Hard spectator scattering 

In this section we consider the spectator scattering mechanism and the calculation of 
{i = 1, 7, 8). The leading corrections from spectator scattering contribute to the branching 
fractions at NLO {0{as)) and are known completely. The NNLO corrections from Qj are 
also known [34], since they can be taken from the heavy-to- light form factor analysis 
in [39,40,42]. In this section we calculate the NNLO corrections from Qg. We find 
agreement with a certain set of logarithmic corrections obtained in [67], and verify the 
important SCET result that contributions at the hard-collinear scale for each operator 
in the effective weak Hamiltonian are taken into account by a universal jet function. To 
complete the NNLO matching calculation for spectator scattering would require results 
for Qi and Q2- This is a rather difficult calculation involving the evaluation of two- loop 
graphs depending on the ratio rric/mb- 

Before presenting our results for Qg, we first review the results for Qy as derived in [34]. 
This will fix some notation and clarify the sub-factorization of t]^ into a convolution of 
hard and jet functions. The calculation makes use of the two-step matching procedure 
outlined in Section [2] to integrate out the perturbative scales ml ^ mfoAqcD- At tree level 
and to leading order in the HQET expansion the result is 

t'l^'\u,u;)= [\rA,C'''^'\T)jf\T,u,u;), (34) 

where 







A,C-.M^=; ,^>^r.u..,.-'-^-l-Hr-u). (35) 

The one-loop correction to the hard-scattering kernel breaks into a sum of corrections to 
the hard coefficient and the jet function according to 

= A7C^^(^) ^ jf + ArC^^W ^ (36) 

where the superscripts denote the (n)-loop correction to each function and the -k denotes 
a convolution over the variable r. Explicit results for each term can be deduced from 
the form-factor analysis in [39, 40, 42] and are listed in the Appendix. Note that while 
the hard coefficient function A^C^^ is particular to the operator Qj, the jet function 
j± is not. It is determined by the matching step SCET/ —>■ SCETn, which contains no 
information about the structure of the operators in the effective weak Hamiltonian at the 
scale rrib. In the SCET description of spectator scattering, therefore, the non-trivial task 
is to determine the corrections at the hard scale m^, contained in the Wilson coefficients 
AjC^^. The contributions at the hard-collinear scale m^A can be obtained by performing 
the convolution in the second term of (IHUjl . 

In what follows we obtain an expression for t^^^^^ in the form flHUl) . derived in the 
following way. We first calculate the hard-scattering kernel directly in QCD factorization. 
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Figure 1: The lowest-order diagram for spectator scattering with Qg. The double-hne 
represents the incoming b quark and the sohd box an insertion of Qg- The photon can be 
attached to any of the four crosses. Only photon emissions from the light quark emerging 
from the Qs insertion contributes at leading power in l/mb- 



but separate the contributions from the hard and hard-coUinear scales using the method 
of regions. We then show that the one-loop contribution from the hard-collinear region is 
exactly AgC^^^^^^j^^. Since both the coefficient function AgC^^^") ~ f /r (with f = 1-r) 
and the jet function j^^^ are non-trivial functions of r, this provides a consistency check 
between the QCD factorization and the SCET formalism, and also a check on our loop 
calculations. The remaining contribution is from the hard region and is identified with 
AgC^^*^^^ Since is a delta function in the variable r, this result is sufficient 

to recover the coefficient function AsC^^'^^\ As mentioned in the Introduction, it is this 
r-dependent function which is needed to obtain the resummed hard-scattering kernel used 
in the numerical analysis in Section [51 



4.1 Qg at tree level 

We start by reviewing the tree-level calculation. The strategy is to evaluate the partonic 
matrix element Ag = {lipi) ^'(^2)7(9) I Qs I Q'ik) b{pb)) and show that it can be written 
in the form i^. The hard-scattering kernel is independent of the exact choice of the 
partonic momenta. We shall work with on-shell quarks in the initial and final states, 
and furthermore work in the reference frame where the perpendicular components of the 
external parton momenta vanish. In this frame, the momenta can be chosen as pi = up, 
P2 = up, k = ujn^/2, ph = rribV = pb — k, and q = E^n+, with u = 1 — u. At leading 
order in 1 / we can write the vector-meson momentum as p ~ mhU- / 2 and the photon 
energy as ^ rnf)/2. The photon's polarization vector lies in the transverse plane and 
is denoted by e±. The power counting is such that uj/rrib ~ A ^ 1. 

The four Feynman diagrams which can contribute at tree-level are represented in 
Figure [TJ The photon can be emitted from any of the four crosses. At leading order 
in A only the diagram where the photon is attached to the light quark produced at the 
flavor- changing weak current contributes. Emissions from the other quark lines are either 
power suppressed or have no projection on the meson LCDAs. For the tree-level scattering 
amplitude one flnds 



(0) _ cQdUs u 1 

nib TT uuu I 
_ nibCF cQdOs u 1 
nib Nc TT uuoj 



u{up) {;^^7.x (1 + 75)} T^u{pb)\ [v{k)llT-v{up) 

^^7.^ ® 7l] ^ t^^7.x ® ll] , (37) 



where u and v represent the free-particle spinor wave-functions and Qd = —1/3 denotes 
the charge of a down-type quark. We have distinguished the MS mass nib{fi) from the 
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pole mass m^, anticipating the one-loop calculation in the next section. To obtain the 
second line we already performed the color trace, so that the notation Fi r2 is to be 
understood as the Dirac structure between quark spinors. In the second line we defined 
the tree-level partonic amplitude A^^\ At tree level there is only one Dirac structure, 
related to the matching of the SCETj operator J^^. 

To proceed further we need definitions for the LCDAs in the low-energy theory SCETn- 
The hght-cone projection operator ^^^(k) onto a heavy state H containing the 6-quark 
is given by 

/I I 
— e"'" (0 [tn_, 0] K^iO)\H{v)) , (38) 

where Qs and are soft and heavy quark fields in HQET; a and (3 are spinor labels. The 
quantity [tri- , 0] denotes a path-ordered exponential along the light cone. Similarly, the 
hght-cone projection operator ^^s{u) onto a light meson state L is defined by 

/r/s 
- e-""+^ (L(p) \^s{sn+) 0] ^,(0) | 0> , (39) 

where the ^ are collinear quark fields in SCET. The hadronic matrix elements of these 
light-cone projection operators, contracted with certain Dirac structures, are the LCDAs 
of the B and V mesons. The exact definitions of the distribution amphtudes needed in 
the analysis are 

J 

(V^(p)|e(.n+) {sn^All^mm = '^^n^V^.{i^ll^\ c^«e— ^0I(k, /.), 

(40) 

where ?7j^ is the polarization vector of the V-meson. 

To extract the hard-scattering kernels we need only the partonic matrix elements. We 
write these as a product of scalar distribution functions multiplied by appropriate Dirac 
spinors. For on-shell matching at leading order in l/m^ the QCD spinors are equal to the 
effective theory ones. At lowest order we have 

$y(a;') = ^''-^^^Up{k)u^{vB-k) = b{u^-J)vf,{k)u^{vB-k), (41) 
$^9'(o) _ 0gr(o) ^^^^p) ^^^^p) _ _ ^^^^p) v^{uv), (42) 

and can be written in the factorized form 



with 



4o) = $''^-'(o)^r8"(°)*$^^"'W, (43) 



^ii(o) / X rribCF eQac^s I , r 

^ 4'^°^{/x7.x(l + 75)},„{7'^-},,. (44) 
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The sub-factorization of tg into the convolution of a hard coefficient with the jet function 
is given by [34] 

AsC '■>{t) = — ; jy{T]u,uj) = -6{t-u). 45 

To show that the hard-scattering kernel 4^ is what appears in the factorization formula 
(in]), we now consider in more detail the hadronic matrix elements of four-quark operators 
in SCETii. Note that the four-quark operator whose hadronic matrix element leads to a 
product of LCDAs has the opposite Fierz ordering ([^^] [gs^t;]) compared to the operator 
whose partonic matrix element matches straightforwardly onto the expression in ( l43l) 
{^hy] [qsi])- In four dimensions the two operators are connected by a Fierz transformation 
according to (see, e.g., [40]) 

Ov^=U^luA'^ + l,)Kqsl''H ^ Oi,^ =e/^^(l + 75)egs^75/i,, (46) 

where the collinear (soft/HQET) fields in each operator are understood to be evaluated 
at different points on the n+(n_) light-cone, and made gauge invariant by inserting ap- 
propriate Wilson lines. In writing ( H6l) we have omitted Dirac structures contributing to 
Oy which have no projection onto the pseudoscalar 5-meson LCDA. Comparing with 
(HOl) . we immediately see that, in the absence of soft-collinear interactions, the hadronic 
matrix element of the operator Oy^ factorizes into a product of (f)\ and (f)^. On the other 
hand, the partonic matrix element of the Fierz-transformed version Oy^ matches fHHj) . up 
to the hard-scattering kernel tg''^^- thus verify that our expression for the tree- level 
amplitude is equivalent to 

We shall perform a similar calculation at one loop in the next subsection. A com- 
plication compared to tree level is that the appearance of IR poles in the dimensionally 
regulated one- loop amplitude prevents one from using the Fierz relation ( l46l) . To extract 
the hard-scattering kernel by comparing the renormalized matrix elements calculated in 
the two Fierz orderings is thus non-trivial, and will require the use of some technical 
details from the SCET analysis used to extract the one-loop jet function in [40,42]. 



4.2 at one loop 

We now turn to a main subject of this paper, the calculation of the one-loop correction 
from Qg to the hard-scattering kernel t^^. The first task is to calculate the amputated 
part of the full set of one- loop Feynman diagrams shown in Figure [2], supplemented by the 
on-shell renormalization factors for the quark fields. Photon emission from the spectator 
quark need not be considered for the case of pseudoscalar 5-meson decay, since the 
four-quark structures appearing in the matching vanish at leading order in 1/mf, after 
projecting onto the meson LCDAs [34]. When a loop integral involves more than one 
scale, we calculate the leading term in the l/rrih expansion using the method of regions. 
All integrals are calculated in dimensional regularization in d = 4 — 2e dimensions and 
with the NDR scheme for 75. The result can be written in the form 

= (/^7.x ® 7l) + 5^'^ (7. ® ll) ■ (47) 

The four-quark structure multiplied by the scalar function A^^^ is related to the matching 
of the operator J^^. It is proportional to the tree amplitude and is used to extract the 
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Figure 2: The one- loop corrections to spectator scattering with Qs- The sohd box denotes 
a Qs insertion and the photon can be attached to any of the crosses. 

hard-scattering kerneL For our choice of external momenta, the function A^^^^ receives 
non-vanishing contributions from the hard, hard-coUinear, and soft momentum regions. 
Since we work with on-shell partonic states, contributions from the collinear and soft- 
coUinear regions are only in scaleless integrals and vanish (this would not be true if 
off-shell regularization were used). We shall label the contributions from the different 
regions as A^i^\A^i^^, and A^P in what follows. Moreover, we define the amplitude A^j^'^ 
to include the as contribution from wave-function renormalization of the 6-quark field, 
which reads 

Zjf -l = -^f- + 31nA + 2^ . (48) 

The renormalization factors for the light-quark fields vanish, to this order in a^. In writing 
the result (H71) . we used the prescription [40] 

Y^l^^^^r^^l^.^lx^7p^ ^ {d-4f{^^Y^^luJ, (49) 

the relevance of which will be explained below. 

The structure multiplied by the scalar function B^^^ is related to the matching of the 
operator J^^. The contributions from individual diagrams contain 1/e poles, but these 
cancel in the sum of all diagrams. Since the matrix element of this four-quark operator 
has no projection onto the 5-meson LCDA, this piece does not contribute to the hard- 
scattering kernel. 

Note that there is no third Dirac structure, which would correspond to a contribution 
from the operator J^. The matching of involves the emission of the n^Ahc component 



15 



of a hard-collinear gluon. Using the equations of motion = i/i-v{up) = ensures 

that the Dirac structure can always be written in terms of the transverse components of 
Dirac matrices. It is easy to see that this precludes taking out a factor of n+Ahc from the 
Wilson line in the operator and attaching it to the spectator quark. 

The function A^^^ has both UV and IR divergences. The UV divergences are removed 
by coupling constant, mass, and operator renormalization (recall that the contribution 
from wave-function renormalization of the 6-quark is included in the definition of A^^^^). 
We define the renormalized parameters as fn^'^^'^ = Zmjnb, etc. As with the vertex term, 
we first compute the QCD amplitude with nf = rih + ni active flavors, and then express 
results in the MS scheme in the n^-flavor theory by renormalizing the coupling as in f|T7j) . 
Doing so, we find that all dependence on rih drops out. The UV renormalized amplitude 
is obtained by making the replacement 

Ad) ^ + aS. = + ( ^^'^+"'(1) + z£) + zg^ - j AW, (50) 

where the various factors at one loop are 

-8^C7^-^J^, Z,,-—-—. (51) 
The M-dependent factor multiplying the counterterm zgj follows from the tree-level coef- 



ficients fl35|45p . To this order in a^, mass renormalization for the 6-quark is needed only 
for the MS mass appearing in the definition of Q^, see (jS]). 

We can extract the one-loop correction to the hard-scattering kernel from the UV 
renormalized partonic amplitude. It is defined by 

^Wio) ^ ^ii(i) ^ ^,,'(0) ^ ^(1) ^ ^(1) _ ^v(i) ^ ^1(0) ^ ^,,-'(0) _ ^Wio) ^ ^ii(o) ^ ^,,'(1)^ (52) 

where the one-loop LCDAs are the renormalized ones. The one-loop contributions to the 
renormalized LCDAs take the form 

= 4^K^^"+4^^^^f, (53) 

^ 4»'*K2' + 4"*K^- (54) 

The renormalization factor Zy^ for the ^-meson LCDA is the Brodsky-Lepage kernel 
[68,69] for a transversely polarized vector meson, and that for the 5-meson was calculated 
in [70]. Here there is an important subtlety, which is discussed in detail in [40,42]. The 
renormalization factors for the B and V meson LCDAs are calculated in the MS scheme 
under the assumption that the light-cone projection operators in ( l38l 1391) are contracted 
with the specific Dirac structures shown in (l40l) . In our case this corresponds to the 
one-loop corrections to the operator O'y^ in fl46p . However, the one- loop amplitude 
is extracted from fl47p and thus multiplies the opposite Fierz ordering, corresponding to 
Ov^ in (H6l) . The Fierz transformation relating the two operators is valid only in c? = 4 
dimensions, and in ci = 4 — 2e dimensions they may differ by terms which vanish in the 
limit d —>■ 4. One can fix these seemingly arbitrary terms by defining a prescription for the 
reduction of evanescent Dirac structures, which ensures that the subtractions in fl3^ are 
performed according to the MS scheme. This is achieved by using the prescription for the 



16 



evanescent operator in (l49l) [40,42]. Since this structure is multiplied by poles related to 
the hard-coUinear region, we can verify that we have performed the correct subtractions by 
checking eq. (l59l) against a corresponding result obtained using the one-loop jet function 
from [40,42]. We describe this cross-check below. 

To evaluate fl32|) it is instructive to rewrite the right-hand side as 



A 



(1) 



+ A 



(1) 



h+hc ' ^-^c.t. 
+^(1) _ ^V(l) 



/,^5'(o)^in(o)^^(i) 



bare 



'k to ^ ^ k: * 



(55) 



We have simplified the above equation by using that the factors Z^^'^ and are delta 
functions, and that the one-loop correction to the bare ^-meson LCDA vanishes for on- 
shell quarks in dimensional regularization (because the integrals are scaleless). Using the 
explicit results from [68-70] to perform the convolutions with the tree-level hard-scattering 
kernel, the first line of (1551) can be written as 



A 



(1) 

h+hc 



+ 4'.' 



(56) 



where 



(1) 




B 




(1) 


Cpa 




Ane 



12 fi 
— + - in — 

e uj 

2\nu] 



5 

2^ 



u 



(57) 



Evaluating (l56l) . one finds that after UV renormalization the IR poles in the sum of 
the hard and hard-collinear regions are exactly subtracted by the poles related to the 
renormalization of the LCDAs. On the other hand, the second line of (1551) vanishes, 
showing that contributions from the soft region are restricted to the S-meson LCDA. 
Therefore, the hard-scattering kernel is free of 1/e poles and insensitive to IR physics. 
This was first verified by the explicit calculations in [67]. 

We have shown that the finite part of the sum of hard and hard-collinear regions is 
t^g^^\ As for Qi, we write the result as 



(1) 



■(0) 



(5^ 



We can identify the finite part of the hard-collinear region with AgC^i^o) * and the 
finite part of the hard region with AgC^^^^'^ ■ For the sum of hard-collinear graphs, 
we find 



he, fin 



a. 



(An) 



[CfJf + CaJa + niif]t 



11(0) 
8 ' 



(59) 
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where n; = 4 is the number of hght flavors, and (Lhc = InmbUj/ fi"^) 

, / 2-21nM , \ 4 
jf = i^hc+ 5 + 21nM Lhc + --12- — 
\ M /MO 

/ 4\ , 21nMlnu , 2 A l\ 4^ . , , 

+ 5 - - lnu + +ln^u 1 + - + -L12 m 

\ u I u \ u I u 



11 , / 1 \ , \ r 76 / 11 1\ , In^M 



f I 1\ fl 1 \ , 2 Inulnw 2 ^ , , 

+ U-^J ^"""U^^J " — — ^ 

2^ 21nM 10 

31 = -Lhe + ^-y. (60) 

Taking the convolution of the one-loop jet function j^^^ listed in the Appendix with the 
leading-order hard coefficient in ( H5l) . we reproduce the above equation. This verifies 
that the sub-factorization of the hard-scattering kernel according to momentum regions 
is equivalent to that in SCET, and also that we have performed the correct subtractions 
in However, we again emphasize that this integrated form cannot be used to obtain 
the resummed hard-scattering kernels used in our numerical analysis in Section [5l 

The finite part of the hard region gives an expression for AgC^^'^'^^ ^ J^'*- In this case 
we have 

AsC^'W ^jio) = ^(1)^ = ^^c^hF + C^M4'^°^ , (61) 



where 



hp = -2L2 -ll + ^-Alnu^L-(l8--] Lqcd + in ( -A + - + 



, ^QCD -r vn - i — 

U J \ U J \ u 

1. 77r2 1^ 27r2 2 4\ 

■15 + — - — + 2 + - + ^ lnM 
12 U 6U \ u u J 

+ (4---^)lnM+(- + ^ + ——] Inulnw- ( 2 + ^ ) In^ u 
u J \u 1 — u) \ J 

+ ( -2 + 7F + 4 + ^in -J Li2(^) + ( - - 4 + ) ^?(«) + ( - - 77^ ) ^(^), 

lu 2[l — u) J \u 2 — uJ \u 1 — uJ 

21nM , \ f I 21nu\ vr^ 3 

21nu — + 21nu L + 4Lqcd + «vr 1 - ^ +2 - — + - 



Inw / , 2 1 \ , / 3 X . , 2 - 

— 1 + — — — mu+ 1 — in u In M — In u 



u \ u j \ 2u 2(2 — u 

+ (4 + I - 2(2^) «(^) + (4 + 2(2^) 

The logarithms L and /^qcd are defined after (12 ip . The terms proportional to L agree 
with the corresponding terms in [67], and can be deduced by convoluting the renormal- 
ization factor z±{t) of the SCET current J^^ obtained in [40,41] with the tree-level hard 
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coefficient AsC^^ in fHSll. We have defined the functions 



In [l-uy{l-y)] 
g{u) = / dy , (63) 

Jo y 

h{u) = /'rfyi^^ll^Ml^. (64) 



These functions have no imaginary part for u G [0, 1] and can be expressed in terms of 
dilogarithms and logarithms, but we shall not give the explicit results here. Since is 
a delta function in r, the result for AgC^^'-^-* is obtained directly from ( l62l) . 



5 Numerical analysis 

In this section we discuss the numerical impact of our results. Our main focus is on 
the branching fractions for B K*'j and Bg 07 decays. The branching fraction for 
B K*'y decays is 



B{B ^ K*7) = [i-^]\A, + A^\\ (65) 



where we have split the contributions from the vertex and hard-spectator corrections 
according to 

A = ^V;;KbX]Ci(/iQCD) A,C^(mb,/iQCD,/i)Ci^i(/^), (66) 

As = ^K:K6X]a(/iQCD)tp(/iQCD,/i)^(^^^0^^//r*0i.l) (/i). (67) 

Results for Bs — > 07 are obtained by making the appropriate replacements. The branching 
fractions depend on a number of parameters, whose values and uncertainties are summa- 
rized in Table [HI The vertex and hard-spectator amplitudes are independently invariant 
under variations of /iqcd and /i. For this reason, their contributions to the amplitudes 
can be studied separately. We discuss each one in turn for the case of i? — > -ftr*7, before 
presenting the final branching fractions for all decay modes in Section 15. 4[ 



5.1 The vertex amplitude 

We begin with the vertex corrections. For these corrections the relevant perturbative 
quantities are the Wilson coefficients Ci in the effective weak Hamiltonian and the SCET 
matching coefficients AjC"^. We calculate the Wilson coefficients in the effective weak 
Hamiltonian using the information summarized in Appendix A of [16] (see also [12, 13]), 
and collect the results for three values of the renormalization scale in Table [21 The 
results for the SCET matching coefficients are simplest when /^qcd = = f^b-, in which 
case all logarithms vanish. We use this choice as our default scheme. Throughout the 
analysis we use the four- loop running coupling with as{mz) = 0.1176, switching from five 
to four active flavors at the matching scale fih = A'-qcd- The vertex amplitude and the 
branching fractions depend rather strongly on (v±- We determine it in Section [5.31 by 
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requiring that the matrix element of Q-j be proportional to the QCD form factor F^^^^ 
at NNLO, finding Cv±{fJ' = ^b) = 0.35 ± 0.05. We will express higher-order corrections to 
the amplitudes in terms of the leading-order result, which is 

= -^V^sVcbC,'' '-^^(ki = -5.48 X 10-^ . (68) 
Up to NNLO, the result obtained using the default parameter values in Table [3] is 

^NNLO 



(0.096 + 0.057i) [a,] + (-0.007 + O.OSOi) [a^] , (69) 



where the first term in parentheses is the NLO (a^) correction and the second term the 
NNLO (a^) correction. The corrections come both from the Wilson coefficients Cj in 
the effective weak Hamiltonian and the SCET coefficients AiC^. Note that we have 
used the effective coefficients in the numerical analysis. This amounts to including 
certain contributions from Q3 . . . Qq, and should be taken into account, if in the future the 
contributions from these operators are worked out systematically. Split into contributions 
from the individual operators, the results for the NLO and NNLO perturbative corrections 
read 

/1NNL0 / \ 

1 = (0.264 + 0.034i) [Qi] - (0.184) [Qr] + (0.016 + 0.023i) [Qg] [as](70) 



+ ( (0.073 + 0.022i) [Qi] - (0.081) [Qr] + (0.002 + 0.008z) [Qs] ] [al 



At both NLO and NNLO the corrections from Qi and Q-j are relatively large, but the 
real parts tend to cancel against each other. Whether this cancellation persists beyond 
the large-/?o limit is an important question. The contribution from to the real part 
of the amplitude is small, but that to the imaginary part is not. It adds together with 
that from Qi to produce a large NNLO correction to the imaginary part. This would be 
a significant effect for CP asymmetries, a topic we leave for future work. 

In Section 15.41 we will study the dependence of the branching fractions on the renor- 
malization scales. To do this we use the SCET Wilson coefficients as given in ( |30ll in 
Section [31 As explained there, this allows us to fix the scale /i = in (v± and study 
the stability of the results under variations in /xqcd and /i/j. Although the expressions 
in the Appendix allow us to vary /iqcd and fih separately, we choose not to do so. For 
simplicity, we set /iqcD = fJ'h and vary them simultaneously. To evaluate the RG expo- 
nents in the SCET evolution factors we distinguish the operators Qj^s and Qi. For Qt.s 
we evaluate the RG exponents using the two-loop anomalous dimensions in a and oj, and 
the three-loop cusp anomalous dimension in the Sudakov factor S. For Qi we evaluate 
the RG exponents using the large-/9o limit. In that case it is consistent to set all SCET 
anomalous dimensions to zero, meaning that we can use the form ( I29l) directly. 



5.2 The hard spectator ampUtude 

The evaluation of the hard spectator amplitude is more complicated than the vertex 
amplitude. It involves a large number of hadronic parameters and the hard-scattering 
kernel contains logarithms of both the hard and hard-coUinear scales. While it is possible 
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Table 2: Wilson coefficients Cj(/i) {i = 1,7,8) at LL, NLL and NNLL. The results at 
NNLL are calculated from the expressions given in [12,13,16], adapted to the operator 
basis in ([21 [3]). The table uses = 4.8 GeV. 





LL 


NLL 


NNLL 


Ci(/i = nib) 


1.11 


1.06 




Ci(/i = y/2mb) 


1.09 


1.04 






1.13 


1.08 




Cf{n = nib) 


-0.312 


-0.303 


-0.294 


Cf (/i = V2mb) 


-0.294 


-0.290 


-0.282 


Cfifi = mb/V2) 


-0.332 


-0.316 


-0.306 


Cf{^ = mb) 


-0.148 


-0.167 




CfiiJ = V2mb) 


-0.141 


-0.159 




Cfifi = mb/V2) 


-0.156 


-0.175 





to fix the scale /iqcd ~ i^b to eliminate some of these logarithms, any choice of the SCET 
factorization scale /i leads to large logarithms in t]^. This can be solved by renormalization- 
group improvement in the effective theory [41]. The hard coefficient AjC^^ is extracted 
at a scale fih ~ /^qcd ~ and and evolved down to the intermediate scale /Xj ~ 1.5 GeV 
by solving the RG equations in the effective theory. The RG-improved hard coefficients 
read [41] 

— e^('^-'^-) / dvU^{u,v,^^H,^^i^)^iC''\v,^^H)■ (71) 

fJ'h/ Jo 

The RG exponents 5* and a are the same as in fl3T]) . The evolution factor U± is the 
solution to the integro-differential equation 

d 

fi—U±{u,v,fih,fJ')= dy'y±{y,u)U±{y,v,fih,fJ'), (72) 
"/^ Jo 

with the initial condition U±{u,v, fih, fJ'h) = 5{u — v). The distribution 7±(?/,m) is the 
anomalous dimension of the operator J^^. A proper treatment of the NNLO matching 
corrections requires this anomalous dimension at two loops, but at present it is known 
only at one loop [40,41]. This adds a small uncertainty to the analysis. The solution to the 
evolution equation is obtained numerically. In the numerical implementation we perform 
the /i-evolution from to /ij in 100 discrete steps. We choose the default renormalization 
scales as /xqcd = fJ'h = fnt and /ij = 1.5 GeV. The dependence on the variable u in 
the resummed AjC^^ is obtained for discretized values of < m < 1. We determine 
the discretization scale by taking more points in u until the numerical convolution of 
the resummed coefficient with the jet function becomes stable. This generally requires 
between one and three-hundred values, although for some cases it is necessary to take 
more values near the endpoints. 

It is natural to evaluate the resummed hard coefficients AiC^^{u, /i) at a scale ~ //j, 
since at that scale the jet function is free of large logs. However, the hadronic parameters 
in Table O are extracted at a low scale = 1 GeV. For a proper treatment one must 
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either run these parameters up to the intermediate scale /Xj ~ 1.5 GeV, or run the hard- 
scattering kernel down to the lower scale. This stage of RG running has been studied 
in [41,70,71]. We have performed this evolution in our numerical analysis but its effect on 
the branching fractions is extremely small. Therefore, in quoting our results, we perform 
the running from fi^ to fii, but ignore that between the scale /ij and the factorization 
scale /i ~ 1 GeV. The shortcoming of this treatment is that the amplitude is not invariant 
under variations of the intermediate scale. However, the dominant effect in this scale 
variation is related the 5-meson distribution amplitude. We account for this in our error 
analysis by assigning a rather large uncertainty to A^. 

A complete treatment of the hard-spectator amplitude is only possible for the NLO cor- 
rections. There are three pieces missing for a full resummed result for the hard-spectator 
term at NNLO: the NNLO hard matching coefficient for Qi, the two-loop anomalous di- 
mension of the current J^^, and the two-loop anomalous dimension of the jet function. 
These missing pieces add uncertainties to the analysis which are difficult to quantify. 
However, we will see that the higher-order corrections from spectator scattering are not 
very important for the branching fractions. 

In addition to the input parameters listed in Table [3], we must also specify the meson 
LCDAs. For the vector mesons we use the Gegenbauer expansion and keep only the first 
two moments: 

(j)y{u) = 6u{l -u)\l + a^{fi)C[^^^\2u - 1) + a^(/x)Cf /^^(2n - 1)1 . (73) 

For the i?-meson LCDA we use the model [72] 

4A^^ u /i 



(cj,/i= IGeV) 



IT UJ^ + fl^ 



/i^ 2{(Tb - 1) , uj' 
in 



(74) 



The i?-meson decay constant in the static limit is 



where to one loop [73] 



KF{fi) = l + — 3 In 2 , V^(/x/„/i) = -— -In (76) 

4:71 \ fi J 2(jo as[Hh) 

We now quote the result for the hard-spectator amplitude to NNLO, accurate within 
the limitations explained above. We find 

/INNLO 

-^f^ = (O.ll + 0.05i) [as] + (0.03 + O.Oli) [a^]. (77) 

Performing the RG evolution of the hard-scattering kernel between fii and the factorization 
scale /i = 1 GeV suppresses the above result by about 10%, or in other words makes 
about a 1% difference on the total amplitude. Split into contributions from the individual 
operators, we have 



^NNLO 



'•hs 



Ah'' 



(0.023 + 0.046z) [Qi] + 0.074 [Qj] + 0.010 [Qs] [««] (78) 



+ (0.004 + 0.003z) [Qi] + 0.025 [Q7] + (0.003 + QMbi) [Qs] [a 
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Unlike the case of the vertex corrections, the individual contributions from the different 
operators are rather small at NLO and especially NNLO. For Qi we have listed the NNLO 
correction found by numerically evaluating AiC^^'^^^ ■ In addition to this correction 
from the jet function, there is also a hard correction AiC'^^'-^-* which is not known. 
Both terms are used for Q7 and Qs- To check the convergence of perturbation theory 
at the intermediate scale Hi ~ 1.5 GeV we split up the contributions from each operator 
into these two contributions. We also separate the NNLO corrections from the Wilson 
coefficients in the effective weak Hamiltonian separate, labeling them with a [w]. For 
these three sources of NNLO corrections, in units of 1/A^^, we have 



Qi : (0.023 + O.Om) [as] + (^(-0.001 - 0.002i) [w] + (0.005 + 0.006i) [jet] J [a^], 
Qr : 0.074 [a,] + 0.002 [w] + 0.015 [jet] + 0.012 [hard]^ [a^], 
Qs : 0.01 [as] + (^Omi [w] + 0.001 [jet] + (0.001 + 0.005i) [hard]^ [a^]. (79) 
In none of the cases is the correction at the jet scale fii = 1.5 GeV unusually large. 

5.3 The SCET soft function 

In this subsection we explain our method for determining the SCET soft function (v±- 
We fix it by requiring that the matrix element of Q7 is proportional to the tensor QCD 
form factor F^^^^ (often referred to as Ti). Using the SCET factorization formula for 
we find 



^0) 



The recent LCSR-based update [25] for the tensor QCD form factor yields F^^^* = 
pB^<j> _ -|- Q_Q4 /iQCD = ^b- Inserting this into (ISUl) and treating the hard- 
spectator term as in the default scheme above leads to Ckl(/^ = i^b) = 0.35 ± 0.05. This 
is considerably smaller than the value — 0.41 used in the SCET analysis in [34], and 
it is mainly for this reason that we find smaller branching fractions below. 

We can use this value for C^y^ to compare the size of higher-order corrections to the 
factorization formula for the form factor. We label the vertex term (v) and the hard spec- 
tator term (hs), and express each as an expansion in a^. Then the individual contributions 
read 

— = (1 - 0.15[a,] - 0.06[a^]) [v] + (0.07[a,] + 0.03[a^]) [hs] . (81) 

For both the as and corrections the vertex term is about twice as large as the hard- 
spectator term and comes with the opposite sign. 

5.4 Branching fractions 

We now convert our results for the amplitudes into estimates for the branching fractions at 
NNLO. The most important uncertainties in the input parameters come from Cv,-, \pz = 
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Table 3: Input parameters used in the calculation of B{B K*'y) and B{B 07). The 
Gegenbauer coefficients in the LCDAs are taken from the LCSR analysis reported in [25]. 



Parameter 


Value 


as{mz) 


0.1176 




-0.040 ± 0.002 


a* JO) 


0.35 ±0.05 


ax(0) 


0.35 ±0.05 


'^bjpole 


(4.80 ±0.10) GeV 


"7.t,pole 


(171 ±2.0) GeV 


/ — / 

^/z = rric/mb 


0.27 ± 0.06 


f 

JB 


/'one; -1- oc;^ tv/TqA/ 
yzyjo It zo J ivie v 


J Bs 


f240 ± 30) MeV 


/(^'Vl GeV) 


fl85 ± 10) MeV 


lore V ) 


(^ 86 ± Q) MeV 


a_|_j^ \L vjreV ) 


U.U4 ± U.Uo 


aYi (l GeV) 


0.0 


af*\l GeV) 


0.15 ±0.10 


a|g(l GeV) 


0.20 ±0.20 


A^^(l GeV) 


(2.15 ±0.50) GeV^i 


aB{l GeV) 


(1.4 ±0.4) 



mc/mf,, Xb, and the renormalization scales. To assess the uncertainty associated with (y^, 
rric and A^, we vary them in the ranges indicated in Table El The scale dependence of the 
branching fraction is completely dominated by the vertex term. We treat this dependence 
as explained in Section ISTTl varying the scale Hh = A^qcd in the range ■mb/\/2 < < 
^J2mh. Including the corrections up to NNLO and discarding terms of O(a^) and higher 
in the branching fractions, we find 

(4.6 ± 1.2 [C^.] ± 0.4 [mj ± 0.2 [A^] ± 0.1 \^\) x 10~^ 
(4.3 ± 1.1 [Ca-] ± 0.4 [m,] ± 0.2 [Ab] ± 0.1 [//]) x 10^^ 
(4.3 ± 1.1 [C<^] ± 0.3 [mj ± 0.3 [Ab] ± 0.1 \^\) x 10"^ (82) 

In cases where the errors are asymmetric, we have taken the average of the higher and 
lower values to get the symmetric form above. The uncertainty in |V^*K6|; which appears 
as an overall factor multiplying the branching fractions, adds about a 10% error to each 
decay mode. To obtain the branching fractions we used the following lifetimes (in units 
ofps) [7] 

r(5°) = 1.527 ±0.008; r(fi+) = 1.643 ± 0.010; r(5,) = 1.451 ± 0.028 . (83) 

In addition to the lifetime differences, our analysis of the three decay modes includes 
differences in the meson decay constants, meson masses, and Gegenbauer moments of the 
light-meson LCDAs (we have assumed that SU(3) violating effects in the i?-meson LCDAs 
are small). Other sources of isospin and SU(3) violation are not included. Concerning 



i3(5+ ^ ir*+7) = 
^(5,^07) = 
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the (f) and K* decay modes, the most important source of SU(3) violation is the difference 
between the SCET soft functions of the two mesons. We discuss this in more detail below, 
giving a result for the ratio of branching fractions of these two decay modes. A study 
of dynamical isospin breaking contributions within QCD factorization was carried out 
in [20]. From this study we expect the dynamical isospin violating effects to make only a 
small difference in the branching fractions. 

It is important to keep in mind that we have not completed the NNLO calculation for 
Qi. The NNLO vertex correction is only an estimate in the large-/?o limit and the NNLO 
hard-spectator correction related to AiC^^ is entirely absent. To study the effects of 
possible deviations from large-/?o limit we assign a 100% uncertainty to the NNLO vertex 
correction from Qi, evaluating the branching fractions using 2AiC^*-^^ and AiC^'-^^ = 0. 
For the hard spectator term we take the NNLO correction as ±1 its NLO value. The 
corresponding uncertainties in the branching fractions, to be added to the errors quoted 
in fl82|) . are ±0.5 for the vertex corrections and ±0.1 for the hard-spectator corrections. 
The uncertainties associated with the unknown corrections to hard spectator scattering 
make little difference for the branching fraction. The uncertainties associated with the 
large-/5o limit in the vertex term are rather large, even though this is an 0{a1) correction. 
We are very conservative with the range in which we vary this correction, but even in the 
only existing calculation of NNLO corrections from Qi beyond the large-/3o limit [53] for 
the inclusive case this is an issue. In that paper the part of the O(a^) correction to the 
matrix element of Qi beyond the large-/9o limit (called P^''^'^'^ [zq) in [53]) remains rather 
uncertain. 

Adding together all the errors mentioned above in quadrature, we obtain the final 
results for the branching fractions 

B{B+^K*+-i) = (4.6± 1.4) X 10~^ 
B{B° K*^-l) = (4.3 ± 1.4) X 10~^ 

B{Bs 07) = (4.3 ± 1.4) X 10"^ (84) 

The NNLO estimates given in (1841) are to be compared with the experimental measure- 
ments summarized in Table 1. We find 

= 1.1 ±0.35 [theory] ± 0.07 [expt.], 
= 1.1 ±0.35 [theory] ±0.06 [expt.], 
= 0.8 ± 0.2 [theory] ± 0.3 [expt.] . (85) 



B{B+ 


-^*^7)sM,NNLO 


B{B+ - 


i^*+7)expt 


B{B^ ^ 


-^*°7)sM,NNLO 


B{BO 


^ A^*07)cxpt 


B{Bs 


07)sM,NNLO 



B{B, ^ 07) 



expt 



Although the results are in reasonable agreement with each other, the theory errors for the 
B — s> A'*7 decay modes are still much larger than the experimental ones. The dominant 
uncertainty is in the SCET soft function C,v^- The remaining uncertainties would be 
greatly reduced by determining the NNLO corrections from Qi to the vertex term beyond 
the large-/9o limit. This would not only directly eliminate the uncertainty in the NNLO 
correction to the hard-scattering kernel, it would also reduce the dependence on the 
charm-quark mass by fixing its perturbative definition. 

Another measurement of interest is the ratio of the branching fractions of the K* and 
decay modes. In the ratio, only the errors in the quantities which are different for 



25 



the Bg, and B, K* mesons add significant uncertainties. Since the spectator scattering 
amphtude is small compared to the vertex term, to a good approximation the error is 
given by that in the ratio Ck*/C<i>- As an example, assuming Ck*/C<I> = 1 ± 0.1, we find for 
the ratio of branching fractions 

—^ = 1.0 ±0.2. 86 

By comparison, the current experimental number is 0.7 ± 0.3. Improved measurements 
of the Bs — > 07 branching fraction, and a more accurate determination of the ratio of 
SCET soft functions, would allow for a comparison between theory and experiment with 
smaller uncertainties than for the branching fractions themselves. 



6 Conclusions 

We computed NNLO corrections to the hard-scattering kernels entering the QCD factor- 
ization formula for B — > decays. We used soft-coUinear effective theory to separate 
contributions between the hard and hard-coUinear scales and to resum large logarithms 
depending on their ratio. For the operators Qj and Qs we obtained exact expressions for 
the hard-scattering kernels for the vertex and hard spectator corrections up to NNLO. 
The results for the vertex corrections provide an explicit demonstration of factorization at 
two loops. For the operator Qi, we estimated its contribution to the vertex correction at 
NNLO using the large- /5o limit. Its complete NNLO correction from hard spectator scat- 
tering was not obtained, but its contribution at the jet scale was evaluated numerically 
and found to be small. 

As an application of our results we provided estimates of the branching fractions for 
B — > K*^ and Bg — > 07 decays at NNLO. The branching fractions are very sensitive to 
the value of the SCET soft function (y^ . We used updated results from QCD sum rules 
for the tensor form factor F^^^^ along with our NNLO results for Qj to find (y^ ~ 0.35. 
Since this value is considerably lower than the default value (y^ ~ 0.41 used in the 
previous SCET analysis in [34], we also find lower branching fractions. Our results for the 
B — > K*^ modes show good agreement with the experimental data, but the theory errors 
are still much larger than the experimental ones. Our result for Bg — > 07, which has a 
comparable theoretical error as in the B K*'~f modes, is also in agreement with the data 
within the large experimental error. The main theoretical uncertainty is in (y^ , which can 
be reduced by improved lattice or QCD sum-rule calculations. On the perturbative side, 
by far the most important issue is the calculation of the NNLO vertex correction for Qi 
beyond the large-/3o limit. This requires the same diagrammatic calculation as the virtual 
corrections to inclusive B Xs'~f, which remains to be done. Our results are also relevant 
for B ^ and B ujj, but for these decays a complete description also requires the 
perturbative corrections to weak annihilation, a topic we leave for future work. 
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7 Appendix 



7.1 Matrix elements 

In the section we give results for the UV renormahzed on-shell matrix elements 

{Q^) ^ {q{p)l{q)mb{p,)) 

in QCD. The results given below are calculated in the MS renormalization scheme with 
Uf — Ufi + rii flavors. For Q7 and Qs we write 

{Qi) — {Q7,tree) 



(87) 

For the results are [43,44] (recall L = In/^/m;,) 

^(1) ^ _^ _ + _2L^-7L- 6.8225 (88) 

-e(1.3333L3 + TL^ + 13.6449L + 13.4779) 

-e^(0.6667L^ + 4.6667L^ + 13.6449L2 + 26.9559L + 26.1412), 



(2) _ I — ■ — I - ■ ■ ------- I 

+5.3333L^ + 36L^ + 96.5798L2 + 144.1712L + 67.6519, 
(2) _ 2.75 3.6667L + 3.5447 4.1546L + 3.4386 



0.5 2L + 2.5 4L2 + 12L + 9.9475 5.3333L=^ + 26L2 + 44.7899L + 27.8816 



-4.8889L3 - 33.9758L2 - 92.3415L - 83.8866, 
(2) _ 0.5 0.6667L + 0.5556 l.llllL + 1.9799 

+0.8889^3 + 6.8889^2 + 19.9050L + 23.8254, 
1.3333L 4L2 + 3.3333L + 0.5483 



D^^^ = + — — + 6.2222L^ + 11.3333L2 + 14.1788L + 0.2934, 

and for Qs we have [45] 

D^s^ = 2.6667L + 1.4734 + 2.0944i + e[2.6667L^ + 2.9468L - 1.1947 + i(4.1888L + 4.1888)] 
+e^[1.7778L^ + 2.9468L^ - 2.3894L - 5.5373 + i(4.1888L2 + 8.3776L + 2.1627)], 

= Di^'^ I -4 - ^^±M ] _ 5.3333L=^ - 32.2802L^ - 50.9612L - 1.8875 



'8F - -8 ^ ^2 

-i(4.1888L2 + 31.4159L + 29.8299), 
D^^l = IS.lllL^ + 31.6617L + 2.38332 + i(23.7365L + 28.0745), 
^ -1.7778^2 -4.0386L- 1.7170 -i(2.7925L + 4.4215), 

= -1.7778L2 -4.0386L + 0.8829 -i2.7925L. (89) 
7.2 The coefficients AiC^^W, AjC^i^i) and j/f ^ 

Here we list the coefficients needed for the numerical analysis of spectator scattering which 
are not written in main text. The lowest order expression for AiC^^ is [34] 

A.c^m^u) = (^) , (90) 
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where 



X 



1 + 4x (^arctanh(Vl — 4x) — i— 



f 1 1 
1 — 4x arctan , , 



(91) 



For the tree-level coefficient was given in (!35l) . The one-loop correction is [39,42] 

vfih vfih vfih 6 



47r 2 



4m 
+ — 

u 



-2 In — - 1 In u + In^ u + Li2(M) 



+2 r^-Tji^ ~T 



4(2 -m) 



u 



-21n-^ - 1 ) Inw + In^M + Li2(u) 

2 



4 4 
+— Li2(M) + - 
uu u 



2 In — - 1 In u + In^ u + Li2(M) 

2 ^ 

' + 41nM-41nM-4 



Li2(M) - y 



(92) 



The one-loop correction to the jet function can be obtained from, e.g., eq. (79) of [40] 
after appropriate replacements. Calling the one-loop correction defined in eq. (79) of [40] 
iY^ after the authors of that paper, we have 

■(1)/ \ ^T^CpCts 1 N 



Nc rrihUJU 



^« -BY/- N 



(93) 



7.3 RG functions 

Here we summarize the perturbative solutions to the RG exponents in ( l30l ITTl) . We define 
the expansion coefficients of the anomalous dimensions and the /3-function as 

2 / (T „ \ 3 



Tcusp \^s) 



^ as ^ /'cts\ -r^ /as\' 



/3(«s) 



-2a, 



^ as ^ { as\^ ^ { as\^ 



47r) 



(94) 



and similarly for the anomalous dimension 7j. In terms of these quantities, the function 
a (and aj with obvious replacements) is given by 



2A 







(95) 



The result for the Sudakov factor S to this same order is 
To f 47r 



4(3^\as{u) 



1 - - - Inr ) + ( ^ - ^ 1 fl - r + Inr) 



/5 



2A 



1 1 2 

— m r 



+ 



47r 





1)'^ 


Uoro 




/?2 






/?oro 



— (1 — r -|- r Inr) + t^tt — — (1 — r) Inr 



To 



(96) 
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where r = as(/i)/as(z/). The cusp anomalous dimension to three loops is 



To 



ye? 









cl 



245 



20 



1347r2 llTT^ 22 ^ , ^ ^ 



27 



45 



418 407r2 



56 



27 



+ CpTpUf I - 



55 



16C; 



27 ^ 



(97) 



and the QCD /3 function is 



11^ 4 



/?2 



{9t 



2857 



Ci + 2C 



205 1415 ^2 

— — Of04 — 
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7.4 Separating scales in the AjC coefficients 

Here we hst the NNLO coefficients AjC"^ in the case where we distinguish Lqcd from L. 
This is achieved by solving the RG equation (l25ll perturbatively, given the form (l22l) for 
the anomalous dimension 7^. This has been done in [74] and we can use those results 
after making appropriate replacements. We find 

A^C^^^) _ ^2 ^2/^4 ^ ^Q^3 ^ 4l2Lqcd + 26.1449^2 + IOLLqcd + 2^2^^ 
+23.5022L + 32.6449Lqcd + 7.8159] 

+CfCa [ - 4.8889L^ - 26.6425^^ - 14.6667LLqcd + 7.33333L^cd 

-63.7859L - 28.5556Lqcd - 83.8866] 

+CFni [0.8889L^ + 5.5556^^ + 2.6667LLqcd - 1.3333L^cd 

+ 17.0161L + 2.8889Lqcd + 23.8254] 

+CFnh (-1.3333L^cD + 2.8889Lqcd - 0.810288) , (99) 



AgC^^^) ^ _^2 [5.3333i^2^Q^j^ + 2.9468L2 + 13.333LLqcd + 16L^cd 
+7.3671L + 43.5941Lqcd + 1-8875 
+2(4.1888^2 + 10.4720L + 20.9440Lqcd + 29.8299)] 

+C^C^ [19.5556LLqcd - 4.4444L^cd + 10.8051L + 20.8566Lqcd + 2.3833 
+i(15.3589L + 8.3776Lqcd + 28.0745)] 

-CfUi [3.5556LLqcd - 1.7778L^cd + 1-9646L + 2.0741Lqcd + 1-7170 
+2(2.7925L + 4.4215)] 

+CFnh [l.7778L^cD - 2.0741Lqcd + 0.8829] , (100) 
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+2/3oLAiC^(^) . (101) 
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